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ABSTRACT 

We report on new sensitive CO J=6-5 line observations of several lumi- 
nous infrared Galaxies (LIRGs: L IR (8-1000/xm) >10 n L Q ), 36% (8/22) of them 
ULIRGs (Lir>10 12 L ), and two powerful local AGN: the optically luminous 
QSO PG 1119+120, and the powerful radio galaxy 3C293 using the James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea in Hawaii. We combine these obser- 
vations with existing low-J CO data and dust emission Spectral Energy Distribu- 
tions (SEDs) in the far-infrared - submillimetre from the literature to constrain 
the properties of the star-forming ISM in these systems. We then build the first 
local CO Spectral Line Energy Distributions (SLEDs) for the global molecular 
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gas reservoirs that reach up to high J-levels. These CO SLEDs are neither bi- 
ased by strong lensing (which affects many of those constructed for high-redshift 
galaxies), nor suffer from undersampling of CO-bright regions (as most current 
high- J CO observations of nearby extended systems do). We find: 1) a significant 
influence of dust optical depths on the high-J CO lines, suppressing the J=6-5 
line emission in some of the most IR-luminous LIRGs, 2) low global CO line 
excitation possible even in vigorously star-forming systems, 3) the first case of a 
shocked-powered high-excitation CO SLED in the radio galaxy 3C 293 where a 
powerful jet-ISM interaction occurs, and 4) unusually highly excitated gas in the 
optically powerful QSO PG 1119+120. In Arp220 and possibly other (U)LIRGs 
very faint CO J=6-5 lines can be attributed to significant dust optical depths 
at short submm wavelengths immersing those lines in a strong dust continuum, 
and also causing the C + line luminosity deficit often observed in such extreme 
starbursts. Re-analysis of the CO line ratios available for submillimeter galax- 
ies (SMGs) suggests that similar dust opacities may be also present in these 
high-redshift starbursts, with genuinely low-excitation of large amounts of SF- 
quiescent gas the only other possibility for their often low CO (high-J) /(low- J) 
line ratios. We then present a statistical method of separating these two almost 
degenerate possibilities, and show that high dust optical depths at submm wave- 
lengths can impede the diagnostic potential of submm/IR lines (e.g. starburst 
versus AGN as gas excitation agents), which is of particular importance for the 
upcoming observations of the Herschel Space Observatory and the era of ALMA. 

Subject headings: galaxies: individual (Arp220) — galaxies: individual (PG 1119+120) 
— galaxies: individual (3C 293) — galaxies: starburst — ISM: molecules 

1. Introduction 

Obtaining unbiased Spectral Energy Distributions (SEDs) for local and distant galaxy 
populations provides a crucial yardstick by which to compare their properties and even- 
tually relate such populations along evolutionary paths within any given galaxy evolution 
framework. Unbiased Spectral Line Energy Distributions (hereafter SLEDs) of the rotational 
transitions of molecules such as CO, HCN, HCO + are of particular importance since: 

• they can trace the mass distribution of the star formation fuel, the molecular gas, 
across its considerable range of properties (n~(10 2 -10 7 ) cm~ 3 , T k ~ (10-200) K) 
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• relative molecular line strengths are in principle extinction-free probes of molecular gas 
properties and AGN-versus-starburst as power sources of IR luminosities of galaxies 
(Meijerink & Spaans 2005; Meijerink, Spaans, & Israel 2006) 

• imaging their emission distribution and velocity fields at mm/submm wavelengths 
yields unique dynamical mass probes of deeply dust-enshrouded star-forming galax- 
ies and QSO host galaxies across the Universe (e.g. Walter et al. 2004; Greve et al. 
2005; Tacconi et al. 2006) 

• cooling and thus the thermodynamic state of the molecular gas (heated by stellar far- 
UV light, cosmic rays, and turbulent motions) is regulated by [CII], [Oil] and high- J 
CO line emission, which in turn may efficiently regulate the local Jeans mass and the 
stellar IMF (Elmegreen, Klessen, & Wilson 2008). 

Much of the star formation in the distant Universe occurs in heavily dust-enshrouded IR- 
luminous systems (e.g. Smail, Ivison, & Blain 1997), challenging to image at optical wave- 
lengths even with the current 8-10 meter-class telescopes. Thus it may well be that in the 
upcoming era of ALMA molecular lines will replace optical and even IR lines as the most 
potent probes of galaxy structure and of the power sources that drive galaxy evolution across 
cosmic epoch. The recent discovery of significant dust optical depths even at submm wave- 
lengths in Arp 220 (a local prototype of dust-enshrouded extreme starbursts) by Sakamoto et 
al. 2008 further highlights how deeply obscured such extreme star-forming systems can be. 

Local templates are key in understanding not only the properties of local galaxy pop- 
ulations, but also of those in the distant Universe. To this end we have undertaken a large 
multi-J CO and HCN line survey of local Luminous Infrared Galaxies (LIRGs) - dust- 
obscured, star-forming systems with SFR~(10-100) M yr _1 (eg. Sanders & Ishida 2004), 
for which the dominant fraction of the bolometric luminosity is in the rest-frame IR with 
Lir(8-1000) /zm)>10 n L . The prodigious star formation events in LIRGs/ULIRGs, most 
often due to dissipative galaxy interactions/mergers of gas-rich progenitors, make them the 
best local analogues of the submillimeter galaxies (SMGs) (Tacconi et al. 2006; Iono et al. 
2009), dust-enshrouded starbursts at high redshifts with even higher star formation rates 
(SFR ~lO 3 M yr~ 1 ) and the sites of a significant part of cosmic star formation history 
(e.g. Hughes et al. 1998; LeFloch et al. 2009). Finally we observed the hosts of two local 
AGN: the powerful radio galaxy 3C 293, and the optically luminous but radio quiet QSO 
PG 1119+120, as part of a pilot study of the molecular gas excitation in the presence of a 
bona fide AGN. The latter remain as the most effective beacons of the most distant galaxies 
where molecular lines have been detected (Walter et al. 2004), and can be the cause of 
distinct molecular gas excitation conditions in high redshift quasars (Schleicher, Spaans, & 
Klessen 2010). 
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In this paper we report on sensitive CO J=6-5 line measurements for galaxies in our 
survey. This allows: a) a first systematic glimpse of their high-excitation molecular gas 
phase, b) constraints on their CO SLEDs, and c) direct comparisons with starburst and 
AGN-dominated systems at high redshifts where predominantly only high-J CO lines are 
currently available (e.g. Solomon & Vanden Bout 2005; Omont 2007). Throughout this 
work we adopt a flat A-dominated cosmology with H = 71 kms -1 Mpc -1 and Q m = 0.27, 
and calculate luminosity distances using the NED calculator developed by Wright (2006). 



2. Sample selection, observations and data reduction 

Our sample of LIRGs where LrR=L(8-1000/im)>10 n L© as estimated from all 4 IRAS 
bands, (e.g. Sanders & Mirabel 1996), of which 36% (8/22) are ULIRGs (L m >10 12 L ), is 
drawn from two CO J=l-0 surveys of LIRGs by Sanders et al. 1991, and Solomon et al. 1997 
(themselves drawn from the IRAS BGS flux-limited sample with f 60/im >5.24 Jy; Soifer et al. 
1987, 1989; Sanders et al. 2003). We then imposed two additional criteria, namely: a) z<0.1 
(the maximum redshift for which the JCMT B-band receivers can tune to CO J=3-2), and 
b) compact CO-emitting regions (sizes obtained from interferometric images in Sanders et 
al. 1988; Planesas et al. 1991; Wang et al. 1991; Downes & Solomon 1998; Bryant & Scoville 
1999; Evans et al. 2000, 2001, 2002) so that one or two telescope pointings can record the 
total line flux up to the highest line frequency of the survey (i.e. 690 GHz of the CO J=6-5 
line). The availability of the CO J=l-0 line flux (Eio/kB~5.5K and n cr it~410 cm -3 ) for all 
selected LIRGs allows an inventory of all the CO-rich molecular gas phases irrespective of 
their average excitation state, and is thus a prerequisite if the fraction of the star-forming 
dense and warm molecular gas is to be determined. The latter is constrained from high-J 
CO line luminosities and reliable global CO (high-J)/ (1-0) line ratios (and hence our second 
criterion for source compactness). 

Finally, along with the aforementioned LIRG sample, two local powerful AGN were also 
included for multi-J CO observations as a pilot study of the possible effects of the AGN on 
the global molecular gas reservoir of their hosts, and to guide future efforts towards a large 
molecular line survey of AGN-selected rather than starburst-selected systems. The two AGN 
chosen are the powerful F-RII radio galaxy 3C 293 whose powerful jets strongly interact with 
the ambient neutral gas (Morganti et al. 2003), and the optically luminous but radio quiet 
QSO PG 1119+120. Both have single dish and interferometric CO J=l-0 line measurements 
that find its emission concentrated in regions with #co^$5"-7" around their active nucleus 
(Evans et al. 1999, 2001, 2005), and thus satisfy our compactness criterion. 
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2.1. Observations and data reduction 

We conducted the CO J=6-5 measurements as part of our multi-J CO and HCN line 
survey of LIRGs in the local Universe. The first CO J=6-5 (691.473 GHz) line measurements 
with the James Clerk Maxwell Telescope (JCMT)Q atop Mauna Kea (Hawaii) were conducted 
for the luminous ULIRG/QSO Mrk231 and the LIRG Arp 193, using the old W/D band 
(620-710 GHz) receiver (operating in SSB mode) on 20 of February and 22 of April 2005 
respectively, under very dry conditions (t 22 oghz ^$0,035). The typical system temperatures 
were T sys ~(3700-5500) K (including atmospheric absorption). The DAS spectrometer was 
used at its widest mode of 1.8 GHz (~780kms~ 1 at 690 GHz), and rapid beam switching at 
^chop=2 Hz with azimuthal throws of 60" resulted in remarkably flat baselines. The beam size 
at 691 GHz is Ohpbw=8". Good pointing with such narrow beams is crucial and was checked 
every 45-60 mins using differential pointing with receiver B3(350 GHz). The latter technique 
can access many more compact sources in the sky in order to conduct pointing checks than 
direct pointing with the W/D receiver, and was found accurate to within cr r ~2.6" (rms). 
The aperture efficiency for these periods was found to be ^*=0.25E 

The main set of CO J=6-5 observations were conducted during several periods in Jan- 
uary, February, March and May 2009 with the upgraded W/D receiver and its new SIS 
mixers (effectively the same type that will equip the ALMA telescopes at this waveband) 
that dramatically enhance its performance. The resulting low receiver temperatures allowed 
very sensitive observations with typical T sys ~(1500-3000) K (including atmospheric trans- 
mission) for T22oghz~0. 035-0. 06. Dual channel operation (two polarization channels aligned 
within <1") further enhanced the W/D band observing capability of the JCMT. The new 
ACSIS spectrometer was used at its widest mode of 1.8 GHz (~780kms~ 1 at 690 GHz), 
and in a few cases (e.g. Arp 220) two separate tunnings per object with an effective band- 
width of ~3.2GHz (~1390kms~ 1 ) ensured adequate coverage of CO lines with FWZI~ 
800-1000 km s -1 . Rapid beam switching at z/ c h op =4Hz (continuum mode) with a throw of 
30" (in azimuth) yielded flat baselines under most circumstances. Finally the pointing was 
checked every 45-60 mins by observing compact sources with W/D receiver itself, and differ- 
ential pointing with the A3 receiver (230 GHz), yielding rms residual error radius of <7 r ~2.2" 
(see Figure 1). The aperture efficiency for these periods is ??*=0.32, and a flux calibration 



1 The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 

2 This is higher than the value adopted by Papadopoulos et al 2007, and we deduced it from a larger 
number of planetary measurements for that period. 
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uncertainty of ~25% (encompassing 77* variations and 3-load calibration uncertainties) was 
estimated from observations of several strong spectral line standards at ~690 GHz. 

A final set of CO J=6-5 observations were obtained recently during several days between 
12 and 23 of January, and 5-6 of February 2010 under dry weather conditions (r 22 oGH Z ^0.05) 
that yielded typical system temperatures of T sys ~(1400-2800) K. The spectrometer setup 
was identical to that used in our 2009 observations with two separate tunnings used to 
cover a wide line in a few objects (e.g. Mrk273, NGC6240). Observations of Mars yielded 
an aperture efficiency of ?7*=0.35, while the overal calibration uncertainty remained within 
~25%. The same beam switching scheme was employed while pointing was checked every 
hour and its accuracy remained within the range shown in Figure 1. 



2.2. Line intensity estimates 

We inspected each 10-min spectrum for baseline ripples and intensity "spikes" in in- 
dividual channels. The edited CO J=6-5 spectra were then co-added and are shown in 
Figure 2, overlaid (mostly) with CO J=3-2 lines. The velocity-integrated line flux densities 
were then estimated from those spectra, after subtraction of linear baselines, using 

Sline = / S,dV = J^K c (x) f TldV = ^^K c (x) f TldV, (1) 

JAY Va, 7rU JAV Va, J AV 

where Fjcmt=15.62 and rf a is the aperture efficiency. The factor K c (x)=x 2 / (1 — e - * 2 ), with 
x — $co/(1-2#hpbw) an d # co =(source diameter) accounts for the geometric coupling of a 
gaussian main beam to a finite-sized disk-like source when 9 co is available. 



2.2.1. Observing compact sources with narrow beams: a bias 

For point-like sources, even with the accurate tracking and pointing achievable by en- 
closed telescopes such as the JCMT, the residual rms pointing errors and the narrow beams 
of large mm/submm telescopes at high frequencies can lead to a substantial systematic re- 
duction of the measured fluxes of compact sources. In Papadopoulos et al. (2009) this is 
shown to be 



GK) = l + 81n2(-^ ) , (2) 
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(see also Condon 2001) where a r /y/2 is the rms pointing error per direction (for identical 
error distributions per coordinate). Thus in Equation 1, when 9 co < 2a T , K c (x) — > G(cr r ), 
since the geometric coupling correction is then overtaken by the pointing error correction. 

2.3. New observations of Arp 220: a revised CO J=6— 5 line flux 

The archetypal starburst Arp 220 has been reobserved during our 2010 observing cam- 
paign, and a significantly larger (~2.6 times) integrated line flux has been measured, along 
with a different line profile. The latter now agrees well with those from high-density gas 
tracers such as HCN J=3-2 (see Figure 3) which places the highest density gas at the high 
velocity component corresponding to emission from the eastern nucleus (see Greve et al. 
2009 and references therein). We thus conclude that the first CO J=6-5 measurement of 
Arp 220 reported by Papadopoulos et al. 2010 has likely been affected by a pointing offset 
of the order of o- r ~4.7" (~3.3" per direction). This is clearly an outlier of the distributions 
shown in Figure 1 but is certainly possible and accounts also for the different CO J=6-5 line 
profile shown in Papadopoulos et al. 2010 (<7 r ~4.7" in this system of two CO-bright nuclei 
~1" apart can yield a maximum apparent emission asymmetry of S(Vi)/S(V2)~0.54 between 
two intrinsically equal emission contributions at the velocity centers Vi and V2 of the two nu- 
clei). Finally our fast-switching scheme allows an independent estimate of the adjacent dust 
continuum from the line-free part of the spectrum, which we find to be S434 /jm =(6.0± 1.5) Jy, 
in good agreement with the value given by Dunne & Eales 2001 rather than that in the 
JCMT calibrator database. The ramifications from this revision are discussed in section 3.2. 



The sources, their coordinates and CO J=6-5 line flux densities (using Equation 1) are 
tabulated in Table 1. The corresponding line luminosities are estimated using 



where ATb is the continuum-subtracted line brightness temperature emerging at the source 
frame, AV and A s are the total linewidth and line-emitting source area, Dl is the lumi- 
nosity distance, and u co ^ cst the line rest-frame frequency. Substituting and converting to 
astrophysical units yields 



2.4. 



CO 3=6-5 line and IR continuum luminosities 




(3) 
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C = 3.25 x 10 7 (1 + z)" 1 ( ^) * ( 2 ( {™ SvdY , ) K km s" 1 pc 2 , (4) 

v 1 V GHz / V m p c / V J y kms ~ / 

while the conversion to ordinary luminosity units (L ), used for the CO SLEDs, is 



L co — 



87rk B z/ 



3 

co, rest 



318 x 104 (iWgik) 1 



lQSKkms-ipc 2 



L - 



(5) 



In Table 2 we tabulate the total IR luminosity L IR of each system, with any contribution from 
a warm AGN-heated dust subtracted when possible. In that table we also list the IR lumi- 
nosity due to star formation l[^ (see section 4), as well as the Rhcn/co = Lhcn(i-o)/^co(i-o) 
and R65/32=Lq Q( - 6 _ 5 - ) /L CO( . 3 _ 2 ^ ) line ratios. The latter is currently the most widely available 
one for conducting comparisons between the excitation conditions in our local LIRG sample 
and high-z starbursts, though for some high-z systems the "neighboring" J=7-6 transition 
rather than J=6-5 is the one available (Tacconi et al. 2006). 



3. Molecular gas in LIRGs: the star forming phase 

The star formation rates of LIRGs are fueled by molecular gas masses of ~ (10 9 -10 10 ) M 
(Tinney et al. 1990; Sanders, Scoville & Soifer 1991; Solomon et al. 1997), at efficiencies 
reaching up to what could be a maximum possible value of SFE=Lt R /M(H 2 )~500 L q /M q , set 
by the radiative feedback of massive stars on the dust mixed with the molecular gas accreted 
by the star forming sites (Scoville 2004). Such high SFEs are one the many indicators 
of the extreme ISM conditions found in such systems, with several studies (e.g. Solomon 
et al. 1992; Aalto et al. 1995; Yao et al. 2003; Greve et al. 2009) revealing also a 
qualitatively different molecular ISM than that of quiescent gas-rich systems. In LIRGs a 
diffuse (n(H 2 ) < 10 3 cm" 3 ) non self-gravitating phase (a likely product of high gas pressures 
inducing an intercloud HI— )>H 2 phase transition, and/or tidal disruption of molecular clouds 
in merging systems) dominates the CO J=l-0, 2-1 line emission. Nevertheless a dense phase 
with n(H 2 )>10 4 cm -3 often contains most of the molecular gas mass present, especially in 
ULIRGs (e.g. Solomon et al. 1992; Gao & Solomon 2004), dominating their CO J+l— >J, 
J+l>3 and HCN line emission. The latter phase is warm, Tkm~(40-100) K, a result of strong 
far-UV photoelectric and/or turbulent gas heating, its high- J CO transitions will thus be very 
luminous and expected to be excellent unobscured "markers" of dust-enshrouded starburst 
regions of galaxies over most cosmic epochs (e.g. Walter & Carilli 2008). 
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3.1. CO 3=6—5 in LIRGs: evidence for dust-affected CO lines 

In Figure 4 the R.65/32 ratio is plotted versus Ljr for all the LIRGs in Table 1, and any 
local or distant galaxies for which this or a similar (high-J) / (low- J) global CO line ratio is 
available. The lightly-shaded area marks the range of R65/32 values expected for the dense 
and warm gas phase that fuels star formation in starbursts, determined using our Large 
Velocity Gradient (LVG) molecular line radiative transfer code for: (n(H2)) = (10 4 -10 6 ) cm -3 
and Tk in =(30-100) K, which encompasses the conditions of star-forming molecular gas. In 
this determination of the R65/32 range we consider only self-gravitating gas, as appropriate 
for any star-forming gas phase, where 



K, 



(dV/dr 
(dV/dr 



obs 



1.54 



virial 



[CO/H 2 

a/«A co 



• (n(H 2 )) 
10 3 cm -3 



-1/2 



(6) 



(e.g. Goldsmith 2001; Greve et al. 2009). Values of K vir ^l correspond to unbound 
gas motions (a~l-2.5 depending on the assumed cloud density profile). The parameter 
A co =[CO/H 2 ]/(dV/dr) along with T^ and (n(H 2 )) classify the LVG solutions (we assume a 
Galactic abundance [CO/H 2 ] = 10~ 4 ). 

From Figure 4 it is clear that for several LIRGs the R65/32 ratio falls well below the range 
expected for star-forming gas, and this happens in some of the most extreme starbursts in the 
local (e.g. Arp220: R 65/32 ~0.2) or the distant (e.g. HR10: R 65/32 ~0.053-0.08) Universal. 
On the other hand all nearby galaxies clearly seggregate into two groups with the star- 
forming ones exhibiting consistently high R65/32 ratios (~0.8-0.9), whilst the SF-quiescent 
ones have much lower CO ratios (~0. 12-0.24) reflecting a much less excited state of their 
molecular gas. Thus in some local LIRGs as well as some of their high-z counterparts, i?65/32 
can be as low or even lower than in SF-quiescent systems. This contradiction between the 
excitation conditions of the star-forming molecular gas and the emergent faint CO J=6-5 
line emission (and thus low R65/32) found in some LIRGs becomes even more acute when 
additional information from other transitions tracing the same phase as CO J=6-5 (e.g. 
those of HCN) is used to further narrow the expected range of R65/32- The (U) LIRGs 
Arp220 and NGC6240 are excellent such examples, with a multi-J line survey of several 
heavy rotor molecules available to constraint the R65/32 ratio (see section 3.2). For all other 
local LIRGs much less such information is available (and is virtually non-existent for their 
high redshift counterparts), with only two small surveys of HCN J=l-0, 3-2 (Gracia-Carpio 



3 If a "neighboring" CO line ratio rather than R65/32 is avalaible for a high-z object, we obtain a range for 
R65/32 by using the available one to obtain LVG solutions over Tk=30-110K, (n(H2))=10 2 -10 7 cm -3 , and 
withK vir ~l. For HR10: R 5 4/2i=0.16^R 65/ 32 -0.053-0.08) 
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et al. 2008; Krips et al. 2008) currently providing the most uniform line datasets probing the 
star-forming phase locally. These yield an average HCN (3-2)/(l-0) brightness temperature 
ratio of (R 32 /io(HCN))~0.55, reaching up to ~0.8-0.9 for the dense gas found in starbursts 
such as M82 (Krips et al. 2008) and Arp220, NGC6240 (Greve et al. 2009). Adopting 
R 32 /io(HCN)=0.55 and K vir (HCN)~l (implemented for an assumed [HCN/H 2 ]=2x 10~ 8 ) 
as constraints on our LVG code yields Tki n =(35-65) K and (n(H 2 )) =(3x 10 4 -10 5 ) cm~ 3 , 
conditions that are indeed typical for the dense star-forming molecular gas (e.g. Guesten et 
al. 1993; Mao et al. 2000) and corresponding to an even narrower range of R 65 / 3 2~0.70-0.85 
(cross-hatched area in Figure 4), excluding several more LIRGs. 

In the absence of dominant low-excitation molecular gas mass in LIRGs (a real possi- 
bility which is discussed in section 5.1), non-negligible dust optical depths at short submm 
wavelengths provide the only other mechanism that can suppress the large CO J=6-5 line 
luminosities expected from their star-forming molecular gas. This has been recently demon- 
strated in the archetypal ULIRG Arp220, where its average dust optical depth of r 8 60/mi~l 
(Sakamoto et al. 2008) has been shown to be responsible for its very faint CO J=6-5 line 
(Papadopoulos et al. 2010). Our current results indicate that such ISM conditions may 
occur often in extreme starburst systems throughout the Universe. 



3.2. The special cases of Arp 220 and NGC 6240 

A recent multi-J HCN, HCO + and CS line survey of these two archetypal ULIRGs by 
Greve et al. 2009 allows placing their CO J=6-5 line luminosities in the best context cur- 
rently possible for such systems. Moreover the significantly higher CO J=6-5 luminosity 
we measured for Arp 220 than that reported by Papadopoulos et al. 2010 necessitates a 
revisit of the analysis for this ULIRG. From the findings of the aforementioned line survey 
it is immediately obvious that for the dense and warm gas phase of (n(H 2 ))~(10 5 -10 6 ) cm -3 
and T k ~(50-120) K dominating the gaseous ISM in these two starbursts: R 65 / 3 2~0.85-0.93. 
This confines it well into the upper (cross-hatched) range of R65/32 values in Figure 4, and 
well above the observed values for both systems. For Arp 220 the new higher value of 
R65/32 =0.21 ± 0.06 remains ~3-6 smaller than that expected for its dense gas phase, and 
thus the evidence for significant dust optical depths at short submm wavelengths quench- 
ing its CO J=6-5 line emission (Papadopoulos et al. 2010) remains strong. Regarding the 
dust emission itself, even for the larger value of S434^ m ~6 Jy and Sg60Atm = (0.55 ± 0.082) Jy 
(Sakamoto et al. 2008), the same analysis as in Papadopoulos et al. 2010 yields T860/im~0.27. 
The latter is smaller than that deduced by Sakamoto et al. 2008, but remains more than 
enough (with a corresponding T434 Atm ~l) to quench the CO J=6-5 line emission in Arp 220 by 
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immersing it into a nearly blackbody continuum at short submm wavelengths. For NGC 6240 
a similar discordance between the R.65/32 expected for its HCN/CS-bright gas and its mea- 
sured R 65 / 32 =0.26 (Table 2) exists. 

Given the large impact that high dust optical depths at short submm wavelengths will 
have to a large number of issues (see discussion in 4.3 and 5), it is worth asking whether a 
large "boost" of the CO J=3-2 rather than a suppresion of the CO J=6-5 line luminosity 
by dust can be responsible for the low R65/32 ratios in some ULIRGs. This can happen only 
if the diffuse non self-gravitating gas phase found in these systems contributes substantially 
to the CO J=3-2 emission (but not to CO J=6-5). For Arp220 the detailed study by Greve 
et. al. yields a diffuse phase with n(H 2 )~(l-3) x 10 2 cm" 3 , T k >40K and A co ~(l-3)xl(r 5 
(kms -1 pc -1 )" 1 , with negligible CO J=3-2 line emission. This is also obvious from the fact 
that even the CO J=2-l transition appears globally subthermally excited in this system. 
Thus a diffuse gas phase cannot be responsible for the very low R65/32 value in Arp220. 
For NGC 6240 this is less clear since its diffuse phase: n(H 2 )~(l-3) xlO 3 cm" 3 , T k >40K 
and Aco~(l~3) x 10~ 6 (kms -1 pc -1 ) -1 can contribute to CO J=3-2 while having negligible 
CO J=6-5 line luminosity. This demonstrates the "degeneracy" between dust-affected and 
genuinely low-excitation CO SLEDs, and the need for multi-J and multi-species molecu- 
lar line observations with excitation characteristics spanning the entire range of molecular 
gas physical conditions in LIRGs. For an earlier such study, conducted for the archetypal 
ULIRG/QSO Mrk231, that finds luminous CO J=6-5, 4-3 line emission from the same gas 
phase as that emitting HCN lines see Papadopoulos, Isaak, & van der Werf 2007. 



4. Global dust emission SEDs in LIRGs: fundamental limitations 

Dust emission SEDs could in principle provide independent constraints on the average 
dust optical depths and thus help settle the issue whether these can affect molecular line 
emission at short submm and far-IR wavelengths. Moreover, particularly well-sampled SEDs 
are available for ULIRGs because of their early importance in deciding their dominant power 
source, i.e. AGN or starburst (e.g. Haas 2001 and references therein). We start by first 
assuming optically thin dust emission SEDs with t\ <C 1 over the entire A~(10-1200) /im 
range, and a superposition of three such SEDs 



1 + z 



Kd (z/)M d , c ( B„(T C ) + ^B„(T SF ) + — — — B^(Tagn) 



M 



d,c 



M 



d,c 



(7) 



This characterizes: a) an AGN- heated dust phase (Tagn, M^agn) (estimated so that its 
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contribution can be reliably subtracted from the global SED), b) a star formation (SF)- 
heated dust (T SF , M d SF ), and c) a cold cirrus-type dust (T c , M d c ). In all cases the important 
850/iin and 1.2mm fluxes have also been included (when available), after been corrected for 
CO J=3-2, 2-1 line contamination. We also correct for any non-thermal radio continuum 
contributions, and a dust emissivity law n^{v) = {v /v Y ^d{^o) with /3=2 is adopted as the 
most suitable for globaldust emission SEDs of LIRGs (Dunne & Eales 2001). The SED fits are 
then used to estimate the total- AGN (L IR ) and SF-related (L^) dust IR luminosities (Table 
2). When a lack of long-wavelength data prevents constraints on the cold dust component, 
the fits are performed by setting T C =22K (the average value obtained from all the other 
SEDs where enough data allowed its estimate). Dust masss and temperatures as well as H 2 
gas mass (obtained from CO J=l-0 data reported in the literature), and the corresponding 
gas/dust ratios are listed in Table 3. For all M(H2) estimates we used a CO-H2 conversion 
factor of X CO =M(H 2 )/L TO(1 _ O) =1M (Kkms^pc 2 )' 1 (Downes & Solomon 1998). 



4.1. Dust in ULIRGs: cold or optically thick? 

From Table 3 it becomes apparent that cold dust mass obtained from typical multi- 
component SED fits always contains most of the mass in LIRGs with M^c/M^sf^-IOO, 
even in extreme starbursts such as Arp220 and Mrk231. This is a direct result of the far- 
IR/submm "excess" found in their SEDs once submm data became available (Lisenfeld et 
al. 2000; Dunne & Eales 2001), and was suspected even when only IR data were available 
because of the high gas/dust ratios deduced with respect to the Milky Way (Devereux 
& Young 1991). Sensitive submm imaging has conclusively demonstrated the presence of 
massive and extended cold dust mass in some LIRGs by spatially separating its emission from 
that of a more luminous SF-heated dust in nuclear regions (e.g. Papadopoulos & Seaquist 
1999, Thomas et al. 2001). From Figure 5 it would then seem that the cold dust mass 
dominates across the entire Lir range of our sample and up to ULIRGs. Nevertheless there 
are some serious problems with this scenario, at least in extreme starburst systems, namely 

• The bulk of the molecular gas mass in most ULIRGs is in a dense, warm, and pre- 
sumably star-forming phase (i.e. n(H2)>10 5 cm -3 , Td us t^40K), and thus most of the 
concomitant dust should also be warm (at such high densities Tk~Td us t)- 

• The M(H 2 )/Mdust ratio obtained from classical dust SED fits of such objects can be 
very low: ~4-20 (Table 3), thus requiring most of the gas in ULIRGs to be atomic 
(with a presumably Galactic gas/dust ratio) so that the global gas/dust ratio can 
attain Galactic values (~140-190, Sodroski et al. 1994). However ULIRGs are found 
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to be exceptionally H2-rich with < 30% of the gas in atomic form (Mirabel & Sanders 
1989, for the X co factor used in our work). 

• Submm imaging of ULIRGs with the SMA (Sakamoto et al. 2008; Wilson et al. 2008) 
found most of the dust emission emanating from very compact regions (<100-150 pc), 
of hot dust (~100-140K), and with M(H 2 )/M dust close to the Galactic value. 

These issues could be resolved if the far-IR/submm "excess" found in the dust SEDs of 
most LIRGs is, in some of them, of different origin rather than cold dust emission. Indeed 
such "excess" can be also generated by emission that is optically thick at far-IR and even 
short submm wavelengths. Recently this has been shown for the compact molecular gas 
disks in Arp220 (e.g. Sakamoto et al. 2008, Matsushita et al. 2009), and is indeed expected 
if most of the molecular gas in ULIRGs is in a very dense state. For (n(Ii2))~(10 5 -10 6 ) cm -3 
(e.g. Solomon, Radford, & Downes 1990; Greve et al. 2009), solar metalicities, and thickness 
of h~(40-60) pc (e.g. Downes & Solomon 1998) such molecular gas disks, if not very clumpy, 
will have dust optical depths of 



r d (A) ~ 1.66 



A 



400 fim 



hMfhOHcosfl)- 1 
10 25 cm -2 



(2.05-30.7) x 



A 



400 /im 
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(6 is the inclination angle, 8=0: face-on). Such extreme gas configurations are the results 
of mergers and/or strong dynamical interactions found in all ULIRGs (Sanders & Ishida 
2004), essentially packing the entire molecular gas supply of two gas- rich spirals within a 
few hundred parsecs. Their emergent IR luminosities will then be less than proportional to 
their dust mass reservoirs, and can be further absorbed by outer dust distributions that are 
optically thick at IR/far-IR wavelengths (Condon et al. 1991; Solomon et al. 1997). For 
such very dense gas disks the dust mass estimates using submm fluxes and temperatures 
obtained from global SED fits (e.g. Dunne at al. 2000) can be significant overestimates since 
an outer cooler dust distribution, while not containing the bulk of the mass, dominates most 
of the observed SED (see also 4.2). This may explain why submm interferometric imaging 
can recover a Galactic value for M(H2)/Md us t in ULIRGs (e.g. Wilson et al. 2008) while 
global SED fits can yield significantly lower ones (Table 3; Dunne & Eales 2001)0. 



4 The M(H2)/Mdust ratios deduced for ULIRGs in Dunne & Eales 2001 must be divided by 5 if an appro- 
priate Xco conversions factor is used for such systems. Then from their Table 10: M(H2)/Mdust(Arp 220)^ 8 
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4.2. Optically thick dust emission: a new degeneracy 



Optically thick dust emission in ULIRGs at far-IR wavelengths has been first noticed 
by Condon et al. 1991 using interfero metric radio continuum imaging to identify the true 
sizes of their starburst regions. To explore its effect on emergent dust emission we follow 
Lisenfeld et al. 2000 and replace the SF-heated and cold dust components inside the brackets 
of Equation 7 with 



where r{u) = t (v / v ) 2 . In Table 3 we list the results of these dust SED fits, and in Figure 6 
show the distribution of the deduced optical depths for all the LIRGs with available submm 
data from the literature. Substantial dust optical depths are found at lOOyum (~4-21) for 
most LIRGs, in accord with earlier results (e.g. Solomon et al. 1997; Lisenfeld et al. 2000). 
Even NGC 1068, a vigorously star-forming LIRG where a massive and extended cold dust 
reservoir withT d ~(10-15) K is revealed via submm imaging (Papadopoulos & Seaquist 1999), 
can have its non-AGN global dust emission fitted equally well with a single-temperature but 
optically thick SED. The latter would then erroneously interpret its cold dust as optically 
thick emission with T 10 o Atm ~4.6. Thus the SED fits using Equations 7 and 9 are degenerate, 
with both cold and optically thick dust emission capable of producing a far-IR/submm 
emission "excess" . This ri R -T d degeneracy can lead to a misinterpretation of the true state 
of the dust in LIRGs, and cannot be broken even for well-sampled SEDs such as those of 
Mrk231 and Arp220 (Figure 7). Only sensitive mm/submm imaging can overcome it by: a) 
spatially separating the cold from the warm and much more luminous SF-heated dust, and b) 
identifying any warm compact dust/gas regions (with potentially significant submm optical 
depths) "nested" inside more extended cooler dust distributions that are optically thick at 
far-IR wavelengths. In Arp 220 such submm imaging has identified hot and spatially extended 
(i.e. starburst-heated) dust reservoirs, with 2r~100pc, T h ~(100-180) K and r h (850;um)~l 
(Eckart & Downes 2007; Sakamoto et al. 2008), within its extended cooler dust distribution 
with T W ~(45-60)K and significant optical depths at the far-IR (r w (100/im)~10; Lisenfeld 
et al. 2000). A simple expression for the emergent SED of such a dust configuration 





1 + z 
4ttD£ 



1 + 



B„(T h ) M d , h F[7hHj 
B„(T w )M d)W F[T w (z/)] e 



k d (i/)M d)W B I ,(T w )F[r w (i/)], (10) 



where F(r) = (1 — e r )/r, t(v)=t (v / V D ) 2 , can easily show the inadequacy of global dust 
emission SEDs in disentangling the true dust components. In Figure 8, after setting T w =40 K, 
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T h =120K r o ,h(850//m)~l, and r OlW (100/xm)~10, M djW =lO 8 M we plot the aforementioned 
SED which makes apparent that hot and compact starburst regions can remain inconspicu- 
ous in the global dust emission SED, apart from an emission contribution in the mm/sub-mm 
domain where a cold dust reservoir would also contribute (and thus the tir-Tj degeneracy). 
This is a result of the high dust temperatures in the compact regions (which makes their 
emission peak at wavelengths where the outer dust distribution remains very absorptive i.e. 
at IR/far-IR wavelengths), and their high dust optical depths over almost the entire observed 
range of a typical SED. The latter diminishes the emission contribution of those compact 
and warm starburst regions well below from being proportional to their dust mass content. 

4.3. The [C II] line luminosity deficit in ULIRGs: a dust optical depth effect? 

One of the most unexpected discoveries made by ISO was the so-called [CII] line lu- 
minosity deficit for some of the most extreme starbursts in the local Universe (Malhotra et 
al. 1997; Luhman et al. 1998, 2003). This was deemed all the more perplexing for ULIRGs 
since the [CII] fine structure line at 157.74 //m is the major ISM cooling line and would 
thus be expected to be particularly bright in such extreme star-forming systems. Various 
explanations have been offered for this such as : a) absorption by a cooler dust foreground 
"screen" , b) saturation of the [CII] line luminosity (while the IR continuum rises unabated) 
predicted for high density PDRs or PDRs in which the local UV radiation field is high and 
the gas density is low, c) soft far-UV radiation fields (see Malhotra et al. 1997; Luhman et 
al. 1998 for a review). Significant sub mm dust optical depths for the bulk of the dust mass 
suggests a simple and common cause for the [CII] and the high-J CO lines luminosity defi- 
ciency in some ULIRGs since in an almost featureless black body dust continuum, no lines 
are expected from a concomitant isothermal gas reservoir. In particular the weakness of the 
CO J=6-5 line argues against the most prominent alternative explanation given for the [CII] 
line luminosity deficit: the dense and/or strongly far-UV illuminated PDRs where in fact 
such high-J CO lines ought to be very luminous. Finally even modest submm dust optical 
depths such as r^o^m^O^ (from Equation 8 and (n(H 2 ))=10 4 cm -3 , h=40pc) correspond to 
T^5g Mm ~1.3 (for (3=2) for the dust continuum at the rest frequency of the [CII] line, yielding 
a F(r) = (l — e~ r )/V=0.56 reduction in line strength with respect to the optically thin case. 

5. CO SLEDs in LIRGs: effects of dust and AGN 

High dust optical depths at short submm wavelengths in some LIRGs (worryingly 
amongst those frequently used as local analogues of high-z starbursts) can significantly sup- 
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press high-J CO or indeed any other high-frequency line emission from their starburst-fueling 
dense molecular gas distributions. In Figure 9 we demonstrate the effects of dust optical 
depth on the CO SLED of Arp 220, the first case where these can be significant (Papadopou- 
los et al. 2010). The highly-excited "intrinsic" CO SLED derived for its massive dense 
and warm gas phase (Greve et al. 2009) is shown overlaid by two SLEDS in which the 
effects of different dust optical depth on the intrinsic CO line luminosities Lj + i j are in- 
serted using Lj^ J =e~ Td ^ J + 1 . J ) x L J+ljJ (valid for an isothermal mix of gas and dust). We 
adopted an emissivity law of (5=2 and r 8 5o Mm =0.2, 1.0 at the rest frequency of CO J=3-2 
(~350 GHz). These values bracket the range of submm optical depths deduced from submm 
interferometry (r 850 ^ m ~l; Sakamoto et al. 2008), and global IR/submm/mm dust SED fits 
(jioo/xm^ 11.6=>- T 8 5o Mm ~0.16; Lisenfeld et al. 2000 for (5=2). 

From Figure 9 is apparent that even the small submm dust optical depths found from 
global dust SED of ULIRGs can have a substantial impact on their CO SLEDs from J=5-4 
and beyond. In such dust-affected CO SLEDs the J-level of the line luminosity peak will 
not reflect a true global gas excitation turnover and thus deeply dust- enshrouded compact 
starbursts can have misleadingly "cool" CO SLEDs. This can particularly affect molecular 
line observations of high-z galaxies where mostly high-J CO lines are currently accessible. 
In such systems, unlike local LIRGs whose molecular gas properties can be constrained by 
many low frequency molecular lines (where dust optical depths remain negligible), separating 
dust optical depth effects from genuinely low-excitation gas can be difficult. 



5.1. Low CO line excitation: the case of Arp 193 

The physical conditions of star-forming molecular gas used to deduce its "intrinsic" CO 
SLED in LIRGs (section 3.1) are typical, however significant deviations do exist. Arp 193 
with a star formation-related IR luminosity of l|| F ' ) =2.2x 10 11 L is notable in having the 
lowest (4-3)/(l-0), (3-2)/(l-0) HCN line ratios: r 43 (HCN)<0.08 (2a), r 32 (HCN)=0.22 (Pa- 
padopoulos 2007; Krips et al. 2008), as well as a faint CO J=6-5 line with R65/32^0.08. These 
are indicative of very low excitation conditions for the bulk of its molecular gas reservoir while 
by comparison NGC 6240, a galaxy with a similar L^^S.Sx 10 11 L , has r 43 (HCN)=0.6 and 
r32(HCN)=0.8 (Greve et al. 2009), indicative of a much denser and warmer gas phase. Radia- 
tive transfer LVG modeling of HCN and CO (6-5)/(3-2) line ratios in Arp 193 (constrained 
also by K vir ~l) are consistent only with low density n(H 2 )~100 cm -3 and warm Tk=(85- 
110) K gas. These are highly atypical conditions for either vigorously star-forming LIRGs 
(where the ISM is usually dominated by dense and warm molecular gas), or quiescent galaxies 
(where molecular gas has similarly low densities but is much colder, Tkm~15 K), and demon- 
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strate that even vigorously star forming systems can have low high- J CO line excitation for 
the bulk of their gas. The corresponding CO SLED for the dense gas in Arp 193 already 
turns over at CO J=3-2, in contrast to the HCN-bright star-forming phase of NGC6240 
which is expected to emit copiously up to CO J=10-9 (Figure 10). It should be noted that 
besides similar lJ| F) , these LIRGs both have M(H 2 )~(0.5-1) x 1O 1O M (for the CO-H 2 con- 
version factor of Downes & Solomon 1998), global SF efficiencies Lj^/M(H 2 )~5OL /M , 
and similar dense gas mass fractions as indicated by Rhcn/co=Lhcn(i-o)/-^co(i-o) (~0.05 
for Arp 193, and ~0.08 for NGC6240). This shows that observables commonly available for 
LIRGs such as IR luminosity, CO J=l-0, and HCN J=l-0 lines, may provide no indication 
of very different molecular gas excitation conditions and CO SLEDs, at least on an individual 
object basis. 

The rapid evolution of starburst events can in principle have short periods during which 
a (U)LIRG stays IR-luminous while its dense gas phase is strongly dispersed and/or depleted, 
by the formation of stars and their concentrated and coherent mechanical feedback (Loenen 
2009). This could place LIRGs such as NGC6240 and Arp 193 at the two ends of such a 
short evolutionary track, but observations of similar objects are needed to decide such issues. 
Finally we note that a population of near-IR selected gas-rich galaxies with ULIRG-levels 
of star formation (~100-150 M yr _1 ) but very low (comparable to the Milky- Way) levels of 
global molecular gas excitation, and star-formation efficiency similar to Arp 193, has been 
recently revealed at high redshifts (Daddi et al. 2008; Dannerbauer et al. 2009). 



5.2. Suppressed high-J CO lines: dust absorption versus low gas excitation 

Once both possibilities of large dust optical depths at short submm wavelengths and 
low gas excitation have been recognized as affecting emergent CO SLEDs, they can be hard 
to tell apart without additional information from other lines tracing the dense star-forming 
phase. These are low-J transitions (thus not affected by dust) of HCN or CS, but can 
be difficult to detect on a routine basis, even in bright nearby LIRGs, since the brightest 
such line: HCN J=l-0, typically has Lg CN ^_ ^ ~( 1/5-1/30) xL' co ^_ ^ (e.g. Gao & Solomon 
2004). An incidious aspect of the aforementioned degeneracy is that some of the most 
extreme starbursts in the Universe can show "cold" CO SLEDs despite hosting intense star 
forming activity simply because their low CO (high-J)/(low-J) line ratios are attributed to 
low-excitation (i.e. non star-forming) molecular gas rather than dust absorption. As already 
discussed in sections 4.1 and 4.2, global dust emission cannot provide a definitive distinction 
either, even with well-sampled dust SEDs in the crucial IR/submm range. 

In Figure 11 we show a possible statistical distinction between these two possibilities 
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using CO J =1-0 (the main bulk molecular gas mass tracer), HCN J=l-0 (the most widely 
used dense gas tracer), along with CO J=6-5, 3-2. The Rhcn/co=Lhcn(i-o)/L'co(i-o) ra tio 
(Table 2) is considered here as an approximate measure of the dense molecular gas mass 
fraction (since n cr it(CO)~410 cm -3 and n cr it(HCN)~2xl0 5 cm~ 3 ), with little sensitivity to 
temperature since for both CO and HCN J=l-0: Eio/ke ~4-5K. The low frequencies of 
~115 GHz (CO J=l-0) and ~87 GHz (HCN J=l-0), ensure that Rhcn/co remains unaffected 
by dust extinction. From Figure 11 becomes apparent that several LIRGs have very high 
R-hcn/co (^ 0.15) and very low R65/32^0.2 ratios, quite unlike typical star-forming galaxies 
with much lower IR luminosities (and thus SFRs). Thus, unless one is willing to entertain 
the possibility of massive amounts of dense but otherwise SF-idle molecular gas in these 
(mostly) merger systems, such LIRGs are likely to be starbursts with dust-suppressed high- J 
CO lines. 

Considering star formation as a process directly fueled by self-gravitating dense gas 
means that star formation switches on rapidly in the dense HCN-bright phase, with its 
typically short dynamical times td yn ~10 5 yrs (much shorter than Giant Molecular Cloud 
evolutionary timescales of ~10 6 -10 7 yrs). This likely yields the tight L ir -L H cn(i-o) correla- 
tion in galaxies (Gao & Solomon 2004), extending "down" to individual GMCs and spanning 
7-8 orders of magnitude in L IR (Wu et al. 2005). It also provides the basis for our aforemen- 
tioned CO-HCN line ratio diagnostic and its ability to "weed-out" dust-affected CO SLEDs 
from genuinely low-excitation ones. In the few cases where multi-J dense gas tracers (e.g. 
HCN, CS) and dust SEDs are available to constrain the state of the dense gas in LIRGs, it 
is always found in the warm and self-gravitating "corner" of the available parameter space, 
typical of star- forming gas (Mao et al. 2000; Papadopoulos et al. 2007; Greve et al. 2009), 
and where R65/32^0.45. Moreover for the lowest possible temperature of Tk~10K in dense 
starless cores deep in GMCs (regulated by cosmic rays rather than photons), a minimum 

R65/32 = E^ U^AbT, _ I ) = 2 (^33.2/T, _ I ) ~ - 32 > (11) 

is expected for the optically thick CO emission (shown also in Figure 4 as a dotted line). Thus 
even for the unlikely case of a high Rhcn/co ratio due to a massive and dense but otherwise 
SF-idle cold gas phase, this is the lowest R65/32 value possible. Values of R65/32^0.3 would 
then imply either a dominant gas phase that necessarily has (n(H2))<10 4 cm~ 3 (and thus also 
low Rhcn/co), o r dust-suppression of the CO J=6-5 line luminosity. In the second case such 
dust-affected LIRGs would then populate the area marked by R65/32^0.3 and Rhcn/co^0.15, 
shown in Figure 11 to be containing several objects. 
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5.3. AGN effects on the high-J CO lines 



The brightest CO J=6-5 lines in our sample are measured in the hosts of two powerful 
AGNs, the optically luminous QSO PG 1119+120 and the FRII radio galaxy 3C293 known 
for a particularly powerful jet (Floyd et al. 2006). These are also the only two objects where 
R 65 /32>1, possible only for CO line emission that is optically thin and remains well-excited 
up to J=6-5, which requires very warm and dense molecular gas. This becomes apparent 
from its maximum value, achieved in the LTE optically thin limit, 



RSS(LTE) = J }" 65, Z k \ (^) = (^) 2 e~ E ^ k = 4e -66.35/T k; (12) 
J(^32,T k ) \ T 32/ V Z/ 32/ 

where J(u, TjJ = (hu/k^) [exp (hzz/keTk) — 1] _1 (the CMB has been omitted for simplicity). 
For PG 1119+120 even R.65/32=2 (i.e. ~(the measured value)-cr) yields T k = 96 K, which is 
rather high for the bulk of its molecular gas. Submm interferometric imaging has recently 
uncovered such high brightness temperatures for the dust (and thus for the concomitant 
molecular gas since T kin >T dust ) in Arp220, but its molecular gas reservoir consists of very 
dense gas with CO line optical depths that are nowhere near the optically thin limit (Greve 
et al. 2009). Comprehensive LVG radiative transfer modeling of such a high CO (6-5) /(3-2) 
ratio, made over a larger gas temperature range of T k =(15-310) K (to allow the possibility of 
very warm gas) yields its best solutions for n(H 2 )~10 5 cm -3 and T k >140K (with solutions 
improving right up to 310 K). Such high gas temperatures could in principle be the result of 
bulk irradiation of the molecular gas by X-rays emanating from a central AGN creating giant 
X-ray Dominated Regions (XDRs), rather than starburst-induced far-UV irradiated Photon 
dominated Regions (PDRs) (see Meijerink & Spaans 2005). Very high CO line excitation 
in AGN hosts that could be attributed to X-ray luminous AGN has been recently observed 
(Papadopoulos et al. 2008), but unfortunately no X-ray data are available for PG 1119+120 
in particular, and the cause of its large CO J=6-5 luminosity remains unknown. 



5.4. CO lines in 3C 293: a shock-powered CO SLED? 

Large amounts of molecular gas (~10 9 -10 10 M Q ) have been discovered in IR-luminous 
powerful radio galaxies locally as well as at high redshifts (e.g. Evans et al. 1999, 2005; De 
Breuck et al. 2005). In 3C293, shocks emanating from a very strong jet-ISM interaction 
(Emonts et al. 2005) may be responsible for a galaxy-wide shock-induced CO line excitation 
(Papadopoulos et al. 2008). A "trademark" of this excitation mechanism could be its 
capability of producing highly excited mid-J and high-J CO lines, even in the absence of 
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high star formation rates, since the shock- induced turbulent heating would deposit most of 
its energy on the gas leaving the concomitant dust mass reservoir much less affected. The 
high luminosity CO J=4-3 and J=6-5 lines in 3C293, with (4-3)/(3-2), (6-5)/(3-2) ratios 
of R43/32=2.3±0.91 and R 65 /32 = l-3±0.54, seem to bear this out by representing some of the 
highest levels of global molecular gas excitation discovered in our sample but in a system 
whose low-level star formation (SFR<4 M yr" 1 ) is insufficient to "power" them. 

We extended our LVG modeling of the extraordinarily high R43/32 and R65/32 ratios in 
3C293 to cover a larger range of Tk=(15-310) K and allow also K vir >l rather than K vir ~l. 
We do so because turbulent heating can drive molecular gas temperatures well above those in 
star forming regions (where typically Tki n < 100 K), while K vir > 1 widens the LVG solution 
search to include non self-gravitating gas phases. These are now more probable given the 
large kinetic energy injection from the powerful jets driving very impressive gas outflows 
in 3C293 (Morganti et el. 2003; Emonts et al. 2005), and the fact that no massive star- 
forming (and thus self-gravitating) gas phase "powering" the luminous high-J CO lines is 
present in 3C 293. The characteristics of the various LVG solution groups are summarized in 
Table 4, where two facts immediately stand out: a) the best solutions are provided by a very 
warm Tk~(130-310) K and dense n(H 2 )~(10 4 -10 5 ) cm -3 gas phase, and b) most solutions 
with Tk<80K are ruled out since they have K vir ^Cl (marked with asterisk in Table 4). 
Such kinematic states are deemed unphysical since gas motions cannot be slower than those 
induced by gas self-gravity (only raising the [CO/H 2 ] abundance by factors of >10 could 
make K vir ~l). On the other hand, almost all of the good solutions found have K vir ~15-49, 
indicating highly unbound gas states. We must also note that in all solutions listed in Table 4 
the fit continues improving well past the Tk=310 K limit of the LVG parameter grid, while 
gas densities remain high (~(10 4 -10 5 )cm~ 3 ). On the other hand the non-AGN cool dust 
SED in 3C293 carries little luminosity (L IR ~3 x 10 10 L o ) and is typical for quiescent ISM. 

Thus in the powerful radio galaxy 3C 293 the molecular gas seems to be in a surprisingly 
hot, dense, and gravitationally unbound state, whilst the bulk of its concomitant dust reser- 
voir maintains lower temperatures and has only a modest IR luminosity. Tellingly, similar 
ISM conditions have been reported for the turbulent molecular clouds in the Galactic Center 
(Rodriguez- Fernandez et al. 2001) where very warm and dense H 2 is concomitant with cold 
dust (Pierce-Price et al. 2000), as well as in 3C326, another another FRII radio galaxy 
(Ogle et a. 2007). In the latter case jet-induced shocks are also suggested as the culprit for 
a LIRG-sized mass (~1O 9 M ) of very warm molecular gas (T k ~125-1000 K), despite star 
formation levels even lower (~0.1M Q yr _1 ) than those found in 3C293. 

In Figure 12 we show the range of CO SLEDs corresponding to the shocked-excited 
ISM environment of 3C293. Clearly more observations of high-J CO lines, now possible 
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with the spaceborne Herschel Space Observatory, are needed to better constrain them and 
verify shock-excited molecular gas over large scales in this powerful radio galaxy and other 
similar systems. The very large molecular gas reservoirs found around radio galaxies at 
high redshifts, and the potential feedback role of their energetic jets in inducing (Klamer 
et al. 2004) or hindering (Ogle et al. 2007) starbursts in galaxy-formation events at early 
cosmic epochs, makes a detailed study of shock-energized molecular gas even more important. 
Distinct and accessible "signatures" of shocked molecular gas states on SLEDs of various 
molecules would be particularly valuable in the upcoming era of ALMA. 



5.5. The diagnostic power of IR/submm ISM lines revisited 

Extremely dense and compact gaseous disks able to attain significant dust optical depths 
at short submm wavelengths in LIRG/AGN systems can have dramatic effects on their 
emergent ISM lines used as AGN-versus-starburst diagnostics. In the presence of such gas 
disks around AGN, ISM line diagnostic using IR and even molecular lines at short submm 
wavelengths (e.g. very high-J CO lines) to discriminate between dust-enshrouded AGN- 
excited XDRs versus starburst-excited PDRs can be become difficult to employ. Relative 
line strengths may then be more indicative of relative dust optical depths rather than intrinsic 
gas excitation properties and their causes. Interestingly all four LIRGs with R65/32 ^ 0.5 
(i.e. IRAS 02483+4302, VIIZw31, PG 1119+120, and Mrk231) have face-on or nearly face- 
on molecular gas disks, indicated either by high-resolution CO interferometry (Downes and 
Solomon 1998) and/or by their narrow CO line widths (3C293 is the only exception). This 
would be expected for gas disks where the smaller dust optical depths in face-on than in highly 
inclined gas disks along the line of sight would enhance their (line)- (continuum) contrast. 
Nevertheless nearly face-on systems with low R65/32 ratios do exist (e.g. IRAS 00057+4021), 
and systematic CO J=6-5 observations of optically-selected QSO samples (which will contain 
many more face-on gas disks) would be valuable in revealing such geometric effects. 



6. Conclusions 

We report on our new sensitive CO J=6-5 observations of LIRGs and two powerful 
AGN with the JCMT, part of a now completed multi-J CO and HCN molecular line survey. 
Our findings are as follows: 

1. Large dust optical depths at short submm wavelengths can suppress the CO J=6-5 line 
in some extreme starbursts by immersing it in strong continuum dust emission. This can 
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yield faint CO J+l— >J lines for J+l>6, even in ULIRGs whose large amounts of dense 
star-forming molecular are intrinsically luminous in such lines. Such high optical depths can 
easily account for the so-called [CII] line luminosity deficity known to exists in such systems. 

2. Similar conditions may be present in high redshift galaxies, yielding deceptively "cool" 
CO SLEDs at high frequencies even in extreme starbursts such as the submillimeter galaxies. 

3. Global dust emission SEDs cannot unambigiously distinguish between high dust opti- 
cal depths at far-IR/submm wavelengths or large amounts of cold dust, as they can both 
contribute to the far-IR/submm part of the dust emission SED in a similar fashion. 

4. The very low CO and HCN line excitation in Arp 193 demonstrates that low global gas 
excitation remains possible even in vigorously star-forming LIRGs. 

5. The suppressing effects of high dust optical depths on the high frequency part of CO 
SLEDs are difficult to distinguish from genuine low gas excitation. Low-frequency (<350 GHz) 
line observations of highly dipolar heavy rotor molecules (e.g. HCN, CS), with their high 
critical densities but unhindered by potentially large dust extinctions at submm wavelengths, 
can "break" this degeneracy. In the simplest such application we propose that LIRGs with 
large HCN/CO J=l-0 but very low CO (6-5)/(3-2) line ratios are likely to have dust-affected 
rather than low-excitation CO SLEDs. 

6. Remarkably high CO line excitation, above that typical for star-forming gas, is found 
for the hosts of two prominent AGN, an optically bright QSO (PG 1119+120) and a radio 
galaxy (3C 293). The latter could be the first known case of a shock-excited CO SLED, likely 
powered by a strong jet-ISM interaction. As a result much of its large molecular gas reservoir 
is hot (Tkin>100K) and dense (n(H2)>10 4 cm -3 ) while most of its dust mass remains "cool" 
(Tdust~15K) typical of quiescent ISM with low star formation rates. 

In summary, the emerging picture of CO SLEDs in LIRGs and AGN seems diverse, with 
high dust optical depths at short submm wavelengths capable of "quenching" high-J CO 
line emission in extreme starbursts while low global gas excitation remaining a possibility 
in such systems. Finally, shocks and possibly X-rays seem capable of surpassing far-UV 
photons as the main excitation contributors in AGN hosts. The now spaceborne Herschel 
Space Observatory is ideally suited for fully characterizing these environments via high-J 
molecular line observations of local IR-bright galaxies and AGN, paving the way for ALMA 
and the study of such phenomena across cosmic epoch. 

We would like to take this opportunity to thank the entire superb crew of people that 
make JCMT a success over so many years. Special thanks to the Telescope System Spe- 
cialists: Benjamin Warrington, Jim Hoge, Jonathan Kemp, and Jan Wouterloot for expertly 
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mission the W/D receiver, as well as Iain Coulson and Antonio Chrysostomou for supporting 
and flexibly allocating the generous amounts of time necessary to conclude this project. We 
would like to also thank Holly Thomas for assisting us with crucial aspects of telescope 
pointing. Like all wines from the island of Santorini, you are all good. Finally we thank the 
referee for useful suggestions that much enhanced this work, and would like to thank again 
the referee of the paper by Papadopoulos et al. 2010 whose suggestions and criticism lead 
to the re-measurement of CO J=6-5 in Arp 220 and the new value of its flux presented here. 
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Fig. 1. — The distributions of pointing offsets obtained during all observing periods and over 
many sky sectors. Offsets with <J<2" are those typical for real source observations (larger 
values are typically obtained after large changes in azimuth and elevation when the telescope 
changed sky sectors between regular observing sessions). During actual source observations 
typically a az ~a ei ~1.6", and a t = + ( x c 2 1 ) 1 / 2 ~2.2". 
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Fig. 2. — CO J=3-2, 6-5 spectra for the galaxies in Table I. The velocity scale is LSR(cz) 
and typical velocity resolutions AV c h~20-35kms _1 . For VIIZw31 the J=2 1 transition is 
also included for more clarity, while for IRAS 08572+3915 the latter transition is the only 
one available. 
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Fig. 2. — (cont.) CO J=3-2, 6-5 spectra for the galaxies in Table 1, with typical velocity 
resolutions AV ch ~20-30kms~ 1 , and in special cases up to 50-60 km s -1 (CO 6-5 in Mrk273, 
CO 2-1 in IRAS 23365+3604). The CO J=6-5 line in IRAS 23365+3604 is shown overlaid 
separately with the J— 2-1 and J=3-2 transitions for clarity, if detected this line seems to 
be present over a narrower FWZI than in the two lower-J transitions. 
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Fig. 3. — The new CO J=6-5 spectrum of Arp220 overlayed on 3xHCN(3-2) from the 
JCMT. Both spectra are clearly rising at higher velocities associated with the eastern nucleus 
that contains the bulk of the high density gas as indicated also by Greve et al. 2009. The 
velocity frame is LSR(cz), and the velocity resolution for both spectra is 28kms~ 1 . 
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Fig. 4. — The CO (6-5)/(3-2) ratio versus non-AGN IR-luminosity (=SF+cold cirrus con- 
tributions, listed as Lir in Table 2) for: a) quiescent (diamonds) and starburst (triangles) 
nearby galaxies, b) our sample of LIRGs (solid boxes), and c) SMGs (stars). We have also 
included a few nearby low luminosity systems (Bradford et al. 2003; Bayet et al. 2006; Mao 
et al. 2000), SMGs at high redshifts (Tacconi et al. 2006; Papadopoulos & Ivison 2002), and 
a Ly-break galaxy (Baker et al. 2004). The star- forming cloud Orion A and the Galactic 
Center (GC, SF-quiescent) are also marked. The shaded areas mark the expected range of 
R-65/32 for typical conditions of the star- forming gas in ULIRGs (cross-hatched) and the less 
extreme star- forming environments of LIRGs (hatched). The dotted line marks the lowest 
possible R65/32 value for the dense and cold GMC cores (n(H 2 )>10 4 cm -3 , and Tk=10K) 
(see 5.2), lower values are possible only for gas with n(H 2 )<10 4 cm~ 3 . 
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Fig. 5. — The cold versus warm dust mass for our sample estimated from typical multi- 
component dust SED fits (Eq. 7). In all LIRGs most of the dust mass is in the cold phase 
Td )C ~(20-25) K, irrespective of IR luminosity, even in extreme starbursts such as Arp220 
and the starburst/QSO system Mrk231. 
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Fig. 6. — The distribution of dust optical depths Ti 00fim deduced by assuming SEDs with 
non-negligible far-IR/submm dust optical depths (Equation 9) versus IR luminosities (with 
the AGN component subtracted when possible: Table 2 second column). The cold dust 
emission and its asscociated masses shown in Fig. 5 are now equally well reproduced by high 
far-IR dust optical depths. 
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Fig. 7. — Top panels: Dust emission of the Arp220 and Mrk231 SEDs using a classic 2- 
component fit (warm dust from SF-sites, and cold cirrus) of its non-AGN emission (section 
4, Equation 7). Bottom panels: Fits using the assumption of significant dust optical depths 
(Equation 9). In the first case large ratios of M diC /M djSF =30(Arp 220),5(Mrk231) are de- 
duced, while for the optically thick fits: Ti O o M m=10(Arp220),4(Mrk231). An AGN hot dust 
emission fit is included in all cases (small curve peaking at short wavelengths). 
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Fig. 8. — The dust SED for a configuration akin to that found in Arp220 (section 4.2, 
Equation 10). The unattanuated hot dust emission from compact starburst regions with 
significant dust optical depths at submm wavelengths (long dashed line), when surrounded 
and absorbed by a more extended and cooler dust distribution with significant optical depths 
at IR wavelengths (its emission indicated by the solid line), remains virtually incospicuous, 
while its peak shifts at longer wavelengths (short-dashed line). The remaining small con- 
tribution of the now attenuated hot component then occurs at long wavelengths where it 
can easily mask as cold dust emission. The emergent total dust SED remains dominated by 
the outer colder dust component (solid line) up to the wavelength beyond which the atten- 
uated emission from the hot dust starts making a significant contribution (marked by the 
continuation of the total SED curve as a dotted line). 
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Fig. 9.— The CO SLED of the ULIRG Arp220, and the effects of dust (see section 5). The 
"intrinsic" (r 850 =0) SLED is produced from the LVG model of the HCN, HCO+ and CS 
ratios reported by Greve et al. 2009. It is obvious that even for the small values of r 850 
deduced for this system (Papadopoulos et al. 2010, this work) the effect on the CO SLED 
beyond J=5-4 can be substantial, making it appear much "cooler" than the intinsic one. 
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Fig. 10. — The CO SLEDs deduced for the HCN-bright (and in principle dense and star- 
forming) molecular gas phase in Arp 193 and NGC 6240, two LIRGs with similar total M(H 2 ), 
Lj^ F \ star-formation efficiency SFR/M(H 2 ), and HCN/CO J=l-0 luminosity ratios (section 
5.1). Surprisingly the physical conditions of the molecular gas in Arp 193 are compatible with 
a total lack of a massive dense (i.e. n(H 2 )>10 4 cm -3 ) gas phase, quite unlike NGC 6240 where 
it dominates its total molecular gas mass and presumably fuels its vigorous star-formation. 
A wider range of CO luminosities for the HCN-bright phase is expected for Arp 193 (denoted 

with thp hart: PYrpnt in trip nh^prvprl CIO transitions! than l\mC! ri94D 
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Fig. 11.— The CO(6-5)/(3-2) (R 65 /3 2 ) versus HCN(1-0)/CO(1-0) (Rhcn/co) brightness 
temperature ratios for our sample (filled circles) and nearby LIRGs (open circles) for which 
these ratios are available. The galaxies with Rhcn/co ^ 0.1 but R65/32 ^$ 0.30 (marked by 
the dotted window) are candidates for having the high frequency part of their CO SLEDs 
suppressed by dust (see section 5.2). 



-41 - 



8xl0 7 



3C293 



6xl0 7 



© 



c 



4xl0 7 



2xl0 7 







high excitation 



lowest excitation 



K . ~1 

vir 



ZL 



cot 


C0( 


C0( 


C0( 


C0( 


C0( 


C0( 


C0( 


C0( 


C0( 


cot 


1—1 

1 


ro 
1 


CO 

I 


i 


cn 

1 


Oi 

1 




□0 

1 


CO 

1 


1— ' 

o 


h- ' 

1 — k 


o 




ro 


CO 




cn 


05 


■v2 


co 


1 

CO 


1 

I— 1 



Fig. 12. — The CO SLEDs possible for the massive and shock-excited molecular gas reservoir 
in the FRII powerful radio galaxy 3C293 where a jet-ISM interaction rather than far-UV 
photons drives the gaseous ISM energetics (see section 5.4). 
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Table 1. LIRGs: Observational parameters and CO J=6-5 line data 



Galaxy a RA (J2000) b Dec (J2000) b D L (z) c / S CO (6-5) dV d 



IRAS 00057+4021 


00 


08 


20 


,58 


+40 


37 


55.5 


194.5(0.0445) 


< 253(1.29,G) 


IRAS 02483+4302 


02 


51 


36 


.01 


+43 


15 


10.8 


225.8(0.0514) 


401 + 131(1. 27,G) 


IRAS 04232+1436 


04 


26 


04. 


.94 


+14 


43 


37.9 


356.4(0.0795) 


< 230(1.21,G) 


IRAS 05083+7936 (VII Zw 031) 


05 


16 


46. 


.51 


+79 


40 


12.5 


239.0(0.0543) 


688 + 207(1. 22,G) 


IRAS 05189-2524 


05 


21 


01 


.41 


-25 


21 


44.9 


186.4(0.0427) 


585 + 200(1. 20,G) 


IRAS08572+3915(NW nucleus) 


09 


00 


25. 


.4 


+39 


03 


54.0 


256.9(0.0582) 


370 + 145(1. 21,G) 


IRAS 09126+4432 (Arp 55 NE) 


09 


15 


55. 


.63 


+44 


19 


58.0 


173.3(0.0398) 


< 500(1.9) c 


IRAS 09320+6134 (UGC 05101) 


09 


35 


51 


,53 


+61 


21 


11.6 


171.1(0.0393) 


< 281(1.35,G) 


IRAS 10173+0828 


10 


20 


00 


.19 


+08 


13 


34.5 


214.4(0.0489) 


< 104(1.27,G) 


IRAS 10565+2448 


10 


59 


18. 


.15 


+24 


32 


34.4 


188.2 (0.0428) 


506 + 148(1. 35,G) 


IRAS 11191+1200 (PG 1119+120) 


11 


21 


47. 


.12 


+11 


44 


18.3 


219.4(0.0500) 


261 + 82(1. 14,K C ) 


IRAS 12112+0305 (NE) f 


12 


13 


46. 


00 


+02 


48 


41.0 


324.4(0.0727) 


327+ 143(1. 22,G) 


LL LL f PITT 7"\ f 

(SW) 1 


12 


13 


45. 


.90 


+02 


48 


39.0 


324.4(0.0727) 


C\C\ / -1 C\ -1 \ 

< 82(1. 21,G) 


IRAS 12540+5708 (Mrk231) g 


12 


56 


14. 


21 


+56 


52 


25.1 


184.1(0.0422) 


1320 + 400(1. 29,G) 


IRAS 13183+3423 (Arp 193) 


13 


20 


35 


,32 


+34 


08 


22.2 


100.2(0.0233) 


< 420 (1.28,G) 


IRAS 13428+5608 (Mrk273) 


13 


44 


42. 


,12 


+55 


53 


13.5 


163.9(0.0378) 


447 + 143(1.45,G) 


IRASF13500+3141 (3C293) 


13 


52 


17 


,77 


+31 


26 


46.1 


194.9(0.0446) 


1086±347(1.29,K C ) 


IRAS 15107+0724 (Zw 049.057) 


15 


13 


13 


,07 


+07 


13 


32.0 


55.0(0.0129) 


904 + 241(1. 14,K C ) 


IRAS 15237+2340 (Arp 220) 


15 


34 


57 


.24 


+23 


30 


11.2 


78.0(0.0182) 


3130 + 810(1.21, G) 


IRAS 16504+0228 (NGC 6240) 


16 


52 


59. 


,05 


+02 


24 


05.8 


104.6(0.0243) 


3321 ±860(1.21, G) 


IRAS 17208-0014 


17 


23 


21 


.92 


-00 


17 


00.7 


186.8(0.0428) 


340 ± 140(1. 17,G) h 


IRAS 23007+0836 (NGC 7469) 


23 


03 


15. 


.60 


+08 


52 


26.3 


69.7(0.0163) 


2355 ± 590 1 


IRAS 23365+3604 


23 


39 


01 


.25 


+36 


21 


08.4 


285.6(0.0644) 


270 ± 95(1. 37,G) h 



a IRAS source name and the most common alternative 
b The source coordinates used for the observations. 

c The luminosity distance in Mpc, and the corresponding redshift in the parenthesis. 

d Velocity-integrated line flux density (in Jykms -1 ) with the (value,type) of any applied corrections 
reported in the parenthesis (see 2.2). All upper limits are at the 2a level. 

c Possible detection for the NE nucleus only, highly uncertain because of a pointing offset of 
~3.85" (thus a beam-shift point-source correction factor of 1.90 has been applied). 

f Tentative detection of the bright NE nucleus (~75% of the CO emission of this double nuclei 
system (Evans et al. 2002). The SW nucleus was not detected, but this particular observation 
may have suffered from a telescope pointing offset. 

g From Papadopoulos et al. 2007 but for n* = 0.25, and a G correction applied (see 2.2.1) 



Table 2. LIRGs: IR luminosities, CO (6-5)/(3-2), HCN/CO J=l-0 line ratios 



Galaxy log (L IR /L ) a log(L};f yL ) a R 65/32 b Rhcn/co c 



IRAS 00057+4021 


11 


.46 


11 


.46 


< 0.24 




IRAS 02483+4302 


11 


.70 


11 


.54 


0.84 + 0.30 




IRAS 04232+1436 


11 


.88 


11 


.87 


< 0.17 




VII Zw 031 


11 


.89 


11 


.69 


0.51 + 0.26 


0.10 + 0.018 


IRAS 05189-2524 


12. 


.06 


11 


.95 


0.57 + 0.20 


0.19 + 0.040 


IRAS 08572+3915 


12. 


.01 


11 


.99 


> 0.6 




Arp 55 (NE nucleus) d 


11 


,27 


10 


.95 


< 0.45 


0.046 ± 0.009 


IRAS 09320+6134 


11 


.88 


11 


.59 


< 0.12 


0.25 + 0.06 


IRAS 10173+0828 


11 


.74 


11 


.62 


< 0.20 




IRAS 10565+2448 


11 


,89 


11 


.81 


0.23 + 0.07 


0.19 + 0.035 


PG 1119+120 


10 


.71 


10 


.46 


3.1 ± 1.2 




IRAS 12112+0305(NE nucleus) 


12. 


,22 


12. 


.08 


0.14 + 0.07 


0.15 + 0.03 


Mrk231 


12. 


.25 


12. 


.14 


0.58 + 0.19 


0.29 + 0.07 


Arp 193 


11 


.57 


11 


.34 


< 0.10 


0.044 + 0.008 


Mrk273 


12. 


.08 


11 


.92 


0.23 + 0.08 


0.32 + 0.05 


3C293 


10 


.49 


10 


.06 


1.31 + 0.54 




Zw 049.057 


11 


,10 


11 


.07 


0.32 + 0.096 




Arp 220 


12. 


.12 


12 


.00 


0.21 + 0.06 


0.18 + 0.03 


NGC6240 


11 


,76 


11 


.58 


0.26 + 0.085 


0.08 + 0.013 


IRAS 17208-0014 


12. 


,32 


12. 


.19 


< 0.10 


0.15 + 0.03 


NGC 7469 


11 


.40 


11 


.21 


0.37 + 0.11 


0.062 + 0.009 


IRAS 23365+3604 


12 


,15 


11 


.79 


0.23 + 0.10 


0.11 + 0.024 



a The SF+cold cirrus (=LiR,(total)-LiR,(AGN)) and SF-related IR luminosities, obtained from 
their respective SED components (see section 4) integrated over A= (8-1000) /im. 

b The R65/32=Lco(6-5)/Lco(3-2) ratio (L' estimated from Equation 3), with CO J=3-2 data 
obtained from Papadopoulos et al. 2010b (all upper limits are at the 2a level). 

c The Rhcn/co=Lhcn(i-o)/^co(i-o) Ta ^° ( c f- Equation 3). The HCN J=l-0 line fluxes are 
from Gao & Solomon 2004 and Gracia-Carpio et al. 2006. In cases of disagreement (Arp 193, 
IRAS 23365+3604, and IRAS 17208-0014) we adopted the values from the latter as more reli- 
able. The CO J=l-0 fluxes used are averages from literature data compiled in Papadopoulos 
et al. 2010b. For NGC 6240 and Arp 220 we used the values reported by Greve et al. 2009. 

d The R65/32 ratio is highly uncertain. We assumed that ~ 2/5 of the total IR luminosity 
of this double-source system emanates from the NE nucleus (same fraction as its CO J=3-2 
luminosity (Leech et al. 2010). 



Table 3. Dust emission SED parameters, gas/dust ratios a 



Galaxy 


Td,sF, T d , c 
(K), (K) 


R b 


M d 
(10 8 M n ) 


M(H 2 ) 
(10 9 M n ) 


M(H 2 )/M d 


r ,T dr (K) 


M dT 
(10 7 M ra ) 


M(H 2 )/M 


IRAS 00057+4021* 


37 


, 6 


134 


25.7 


4.1 


1.5 


8.7, 65 


0.4 


1057 


IRAS 02483+4302 


36, 


20 


17 


4.6 


3.6 


8 


11, 49 


2.0 


176 


IRAS 04232+1436 


35 


, 9 


33 


23.7 


10 


4 


24, 69 


0.80 


1220 


VII Zw 031 


51, 


23 


80 


3.6 


11.5 


32 


21, 55 


11 


107 


IR 05189-2524 


60, 


28 


33 


1.1 


3.9 


35 


4.8, 72 


0.7 


557 


IRAS 08572+3915 


52, 


19 


21 


8.2 


1.6 


1.9 


10, 85 


6.5 


24 


IRAS 09320+6134 


40, 


25 


19 


2.9 


4.8 


16 


5, 47 


3.7 


131 


IRAS 10173+0828 


48, 


23 


27 


1.5 


6.0 


39 


10, 64 


0.60 


993 £ 


IRAS 10565+2448 


41, 


23 


7 


1.9 


6.4 


32 


6.2, 61 


1.3 


516 


PG 1119+120 


41, 


22 


29 


0.3 


0.5 


16 








IRAS 12112+0305 (NE nucleus) 


43, 


25 


10 


3.3 


7.5 


22 


5.3, 60 


2.8 


272 


Mrk231 


44, 


28 


5 


1.9 


7.0 


36 


4, 66 


2.1 


336 


Arp 193 


42, 


27 


10 


0.7 


4.7 


59 


0.5, 38 


3.9 


119 


Mrk273 


52, 


29 


18 


1.3 


5 


38 


4, 65 


1.3 


385 


3C293 


57, 


22 


459 


0.2 


4.6 


177 


5, 38 


0.3 


1398 


Zw 049.057 


33, 


17 


4 


0.6 


0.9 


13 


7, 52 


0.4 


214 


Arp 220 


46, 


22 


30 


5.7 


6.0 


11 


10.5, 62 


1.8 


339 


NGC6240 


52, 


26 


37 


8.0 


8 


10 


5.6, 59 


0.8 


976 


IRAS 17208-0014 


46, 


24 


19 


5.3 


13 


24 


7.4, 60 


3.2 


416 


NGC7469 


42, 


27 


9 


0.4 


3.5 


72 


4.5, 55 


0.6 


602 


IRAS 23365+3604 


60, 


30 


84 


8.2 


7 


8 


3, 54 


3.0 


223 



a The full dataset used in these SED fits will be presented in an upcoming paper (Papadopoulos et al. 2010, in preperation). 

b The Md,c/M di sF cold/warm dust mass ratio. 

*Very sparsely sampled dust SED, highly uncertain fit. 
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Table 4. Physical conditions of the molecular gas in 3C 293 



(4-3)/(3-2), (6-5)/(3-2) a 


T k 


n(H 2 ) 


K • b 

J-V vlr 


x 2c 




(K) 


(cm' 3 ) 






2.3, 1.3 (measured) 


15-125 


(1-3) xlO 5 


0.027-0.154* 


2.21-1.6 




130-310 


10 5 


15.4 


1.6-1.23 


1.27, 0.7 (minimum) 


15-30 


(0.3-3) xlO 5 


0.027-0.154* 


1.16-0.74 




40-55 


3 x 10 3 


0.27* 


0.7 




75-115 


3 x 10 4 


27.4 


0.63 




135-255 


10 4 


4.9-15.4 


0.68-0.59 




260-310 


10 4 


49 


0.58-0.51 


2.3? 2.34 (maximum) 


15-125 


(1-3) xlO 5 


0.027-0.154* 


2.42-1.94 




130-310 


10 5 


15.4 


1.93-1.54 


1.27, 1.3 (medium) 


15-125 


(1-3) xlO 5 


0.027-0.154* 


1.75-0.89 




130-310 


10 5 


15.4 


0.88-0.21 


1.27, 2.3 (medium) 


15-125 


(1-3) xlO 5 


0.027-0.154* 


2.11-1.56 




130-310 


10 5 


15.4 


1.56-1.15 


2.3 f 0.7 (medium) 


15-30 


(0.3-3) xlO 5 


0.027-0.154* 


1.66-1.4 




40-55 


3 x 10 3 


0.27* 


1.36 




80-130 


3 x 10 4 


27.4 


1.31-1.33 




145-245 


10 4 


4.9-15.4 


1.35-1.30 




250-310 


10 4 


49 


1.30-1.25 



a CO line luminosity ratios (L^ ± a x )/(h' y ± a y ) (Eq. 3), with (type) denoting 
the L ± a values used, for example: (measured) = (L 43 /L 32 , L 65 /L 32 ), 
(minimum)= [(L^ - t74 3 )/( L 32 + a 32 ), (Lg 5 - a 65 )/ (L 32 + a 32 )] etc 

b The minimum virial parameter of the LVG solutions (Eq. 6 for a — 1). 
Values K vir < 1 are forbidden (see 3.1) and are marked with an asterisk. 

c The x 2 = Ek VK) 2 [Rk(LVG) - R k (obs)] 2 of the LVG fit of the observed 
line ratios Rk(obs), when a range of x 2 values is reported, the low values always 
correspond to the highest temperatures. 

d For R43/32, the observed value is so high that is effectively identical to the 
LTE, optically thin limit: (4/3) 2 = 1.77, thus we do not consider the +a case. 

*Kinematically unphysical solutions (gas motions slower than virial). 



